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We suggest that the phase transition observed in NaV205 at T = 34^" is not a spin-Peierls transition, 
but a charge ordering transition, related to the formal presence in this system of equal number of 
V*+ and ions. Below Tc, ions form a zigzag structure, which is consistent with the 

experimentally observed doubling of the lattice period in a and h directions. We show that this 
charge ordering also results in the alternation of spin exchange constants along the ^-direction, 
which opens a gap in the spin excitation spectrum. We emphasize the role of lattice distortions 
around V ions both in the formation of the charged ordered state and in the spin-gap opening. 

PACS numbers; 64.70Kb, 75.10Jm 

In this paper we discuss the interplay between charge 
and spin degrees of freedom in NaV205, which initially 
was identified as an inorganic spin-Peierls (SP) material 
similar to CuGeOs. Indeed, according to the first X- 
ray studies of this compound Q V'*^ ions form spin-^ 
chains, which are separated from each other by nonmag- 
netic V^^-chains. Below Tc = 34K, magnetic susceptibil- 
ity. X-ray and neutron scaterring measurements indicate 
the opening of a spin gap accompanied by doubling of 
the lattice period in the chain direction |^^. 

However, the situation in NaV205 is evidently more 
complicated and more interesting than in conventional 
SP systems. First of all, the original crystal structure 
was recently questioned l^-Q]. According to the new X- 
ray and neutron data, the structure of NaV2 05 was iden- 
tified as a centrocymmetric group Pmmn, in contrast to 
the earlier assignement P2imn Q. The new structure 
implies that all V sites are equivalent. As an average 
valence of V in NaV205 is 4.5-I-, there is one d-electron 
per two equivalent V ions, which makes the description 
of this compound as a spin-^ chain material not obvious. 

Other indications that the physics of NaV205 may be 
different from that of, e.g., CuGeOs come from thermo- 
dynamic data. In particular, the study of specific heat 
Igl shows that the entropy of the transition in NaV205 is 
larger than the entropy of a pure SP system. The ratio 
where Aq is the value of the spin gap at T = 0, is for 
NaV2 05 '--^ 6 - much larger than for the known SP mate- 
rials, where it is close to 3.5, which follows from the weak- 
coupling mean field theory of SP transition [|o| . Another 
evidence that the phase transition in NaV205 may be of 
different nature comes from a rather weak dependence 
of Tc on magnetic field: the shift of Tc is ~ 5 times 
smaller [^,|| than the theoretical predictions [ p^HTl| ] . Yet 
another spectacular difference in the behavior of NaV205 
and CuGeOa is the temperature-dependence of the ther- 
mal conductivity jl^ : below Tc the thermal conductivity 
of NaV205 increases by a factor of 5, while in CuGcOa 
the corresponding anomaly is very weak. 
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FIG. 1. Schematic representation of the crystal structure of 
V-O plane of NaV2 05. Oxygenes are located at the corners of pla- 
quettes and vanadium ions are located at their centers; the shaded 
plaquettes are vacant. Also shown are the relevant d^iy-orbitals of 

V ions and vanadum ladders (dashed lines). 

All these facts show that the phase transition in 
NaV205 is not an ordinary SP transition. We suggest be- 
low that the main phenomenon responsible for the tran- 
sition is, in fact, a charge ordering (CO), related to the 
formal presence in this system of equal number of 
and ions. While at high temperature there is a rapid 
oscillation between the V^^ and states, so that on 
the time scale larger than the electron hopping time all 

V ions look equivalent {V^'^^), at low temperature a CO 
transition, similar to the Verwey transition in Fe304 [ p^ , 
occurs. The opening of a spin gap, acompanying this 
transition, is just one of the consequences of the charge 
ordering. The appearance below Tc of two ineqivalent 

V sites, which were crudely identified as and V^+, 
was indeed observed in the recent NMR experiment [Q . 
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However, the detailed nature of this CO and its relation 
to the opening of the spin gap remains still unclear. Elu- 
cidation of this relation is the main subject of this paper. 

The structure of NaV205 consists of V2O5 layers with 
approximately square oxygen lattice. Two third of the 
oxygen plaquettes are occupied by V ions located in the 
middle of the plaquettes (see Fig. 1). Due to details of the 
crystal structure, only one of the three t2g-orbitals of V, 
namely the d2;y-orbital, is occupied |^. As the d-d hop- 
ping goes mainly via intermediate oxygens, the diagonal 
(or next-nearest-neighbor (nnn)) hopping amplitudes, tj_ 
and are much larger than the nearest-neighbor (nn) 
hopping amplitude t^y, which makes NaV205 a quater- 
fiUed system of vanadium two-leg ladders (See Fig. |^). 
Tight binding fitting of the LDA band structure gives 
tjL ~ 0.35eV, t|| - O.lSeV, and tj:y - 0.012eV g, aU 
remaining hopping amplitudes being much smaller. Us- 
ing the empirical Harrison rules Horsh and Mack ob- 
tained similar values for t±_ and but a much larger 
~0.3eV [0. 

With the parameters of Ref. 1^ one can understand the 
one-dimensional spin structure of NaV2 05. Each rung is 
occupied by a single electron in a bonding state with the 
energy —t±. For strong on-site Coulomb repulsion the 
energy cost of electron transfer on a neighboring rung 
is (at least) 2t±. For 2t± > At\\ the electron hopping 
between rungs is energetically unfavorable and can oc- 
cur only virtually. Thus, although the charge degrees of 
freedom are not frozen and all V ions are equivalent, the 
system remains insulating. In Refs. NaV205 was 

described as an effective spin system with spins being 
localized on the rungs of vanadium ladders (instead of 
vanadium sites) . As the spin-exchange between electrons 
of neighboring ladders is weak Pjl5|], the spin system 
is quasi-one-dimensional (each ladder corresponding to a 
spin chain). 

In this paper we go beyond the decription of NaV2 05 
as a spin-chain material and include also the charge dy- 
namics. We start from the electronic Hamiltonian, which 
includes the electron hopping in ladders (with the ampli- 
tudes t± and <n) and between ladders (with the amplitude 
txy) [see Fig. as well as the Coulomb interaction be- 
tween electrons on different V sites: the on-site Coulomb 
repulsion U is the largest parameter and is here taken in- 
finite, the nn interaction V2 (e.g., between sites C and D 
on neighboring ladders [see Fig. ^), the nnn interaction 
Vi inside ladders (e.g., between sites A and E or B and 
E [see Fig. H), etc. 

Next we make a projection on the subspace of states 
with one electron per rung. Four different states of a 
single electron on a rung can be represented as the eigen- 
states of spin 5* = ^ and isospin T = ^ operators, 
\^S^) (8) I^T^), where = ±i corresponds to an up- 
per/lower position of the electron on a rung (see Fig. |^). 

The effective spin-isospin Hamiltonian can then be 
written in the form: 



H = Ho + Hi. 
Here, Hq is a pure isospin Hamiltonian: 

r r ^ ^ 



(1) 



(2) 



where the vector r runs over the sites of an effective lat- 
tice located on the rungs of V ladders. The first term 
in Eq.(H) describes the electron hopping in the rungs, 
the second term is the nnn Coulomb interaction between 
electrons occupying neighboring rungs r and r -I- b of a 
ladder (b is a unit lattice vector in the ladder direction), 
the third term is the repulsion between nn electrons from 
neighboring ladders and Hq stands for all other interac- 
tion terms (the longer range Coulomb terms, the inter- 
action via the lattice, etc.). 

The Hamiltonian Hi in Eq.(^ decsribes the interac- 
tion between spin and isospin degrees of freedom. The 
strongest spin-isospin interaction occurs inside ladders. 
The corresponding Hamiltonian for a single ladder has 
the form: 



HsT — — 



4f2 



- - S„- S„+iY- + T„- T„+i - 2T.^T^j^-^ 
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— + S„-S„+i II — — T„-T„_|_i 



(3) 



where n runs over rungs of the ladder. 

The Hamiltonian Eq.(|) originates from the virtual 
electron hopping between two neighboring rungs of a lad- 
der and A is the increase of energy due to electron trans- 
fer on the neighboring rung. The value of A is deter- 
mined by the isospin Hamiltonian Hq and depends on a 
state considered. If, for instance, we neglect the Coulomb 
terms in the disordered phase, then A = 2tj^, while for 
the extreme zigzag CO discussed below: A ~ 2Vi. Note, 
that in the U ^ 00 limit A is finite. Since the hopping 
between ladders, t^y, is much smaller than t± and i[| [^, 
we ingnored the corresponding contributions to Hi . 

The spin-isospin structure of the exchange interaction 
is a consequence of the Pauli principle: the wave function 
of two neighboring electrons in the intermediate state 
(when they occupy the same rung) should be antisym- 
metric under exchange of both their spin and isospin co- 
ordinates. Therefore, if the total spin of two electrons on 
neighboring rungs, S = S„ + S„+i, is zero (i.e., the spin 
state is antisymmetric), the isospin state of these elec- 
trons has to be symmetric: T — 1, where T = T„-f T„+i, 
and, vice versa, if = 1, then T = 0. Correspondingly, 
the first term in square brackets of Eq. ||is the projection 
operator on S' = 0, T = 1, = state, while the second 
term is the projection operator on 1, T = state. 
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The projection of the total isospin of two electrons, T^, 
in both terms is zero, because we assumed the on-site 
Coulomb repulsion U to be infintely large, so that the 
hopping between rungs is only possible if electrons are 
located on different chains of the ladder. 
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FIG. 2. Two types of CO: the zigzag structure (a) and the old 
structure, in which chains are separated by V^"*" chains (b). 

Large circles denote ions with a large ionic radius, while small 
circles denote small nonmagnetic ions. Arrows denote shifts 

of oxygen atoms around ions in one ladder (for simplicity, the 
shifts induced by V'*"'' ions in other ladders are not shown). 

Since TJi, resulting from the virtual hopping between 
rungs, is smaller than i/oi we can start by considering 



only the charge degrees of freedom first. The isospin 
Hamiltonian describes the competition between the 
hopping along rungs, which tends to make V ions equiva- 
lent, and the Coulomb interaction, which favors a CO. If, 
for instance, we leave only the first two terms in Eq.(^, 
then the ladders become decoupled and the Hamiltonian 
of each ladder is the Hamiltonian of the Ising model in 
perpendicular magnetic field. For V\ < 2ti_ its ground 
state is disordered (all vanadium ions are equivalent), 
while in the opposite case, isospins in the ground state 
are ordered antiferromagnetically: (2T^) = r7(— 1)". This 
CO corresponds to the zigzag occupation of vanadium 
sites by electrons as shown in Fig. ||a , where the large 
circles denote V^^ ions with a large ionic radius, while 
small circles denote small nonmagnetic ions. (Of 
course, if the ladders would be really decoupled, the sign 
of the order parameter in each ladder is arbitrary.) 

On the other hand, the third term in Eq.(^) favors 
the structure shown in Fig. [^b , which would give one- 
dimensional spin chains of V^^ {d}, S = ions and is 
the initially assumed crystal structure If we would 
take into account only the second and the third terms 
with V2 = V2Vi, the zigzag structure (Fig. ^a) would 
have lower energy than the chain structure (Fig. ^b). Of 
course, it is necessary to include all longer range Coulomb 
interactions to find the real ground state of the system. 
The calculation of the Madelung energies of these two- 
dimensional structures shows, that they are rather close, 
the chain structure of Fig. ^ being slightly more favor- 
able. 

There is, however, another type of interaction between 
charges - the interaction via lattice distortions, which, 
unlike the Coulomb interaction, is local. As shown in 
Fig. ||a, this latter interaction definitely favors the zigzag 
structure: large V^~^ ions push out neighboring oxygens 
(the directions of oxygen displacements in a ladder are 
shown in Fig. ^ by arrows). As a result, neighbor- 
ing plaquettes in a ladder would preferably be occupied 
by smaller V^^ ions. Roughly speaking, the lattice re- 
laxation around empty and occupied vanadium ions en- 
hances the nnn repulsion, which favors the zigzag order- 
ing. 

Note, that the zigzag structure shown in Fig. ||a im- 
mediately gives us a new periodicity with the doubling 
of periods along both a and b directions, in accordance 
with the experiment Thus, one does not need to in- 
voke the SP mechanism to explain the distortion pattern 
observed below Tc. We also note, that there exist four 
equivalent realizations of the zigzag structure of Fig. |^a, 
which differ by the location and orientation of pairs of 
nn V*^ ions (large circles in Fig. |^a). This may provide 
a clue to understanding of the four-fold increase of the 
period in c-direction below . 

The smallness of Tc = 34if in NaV205 as compared to 
the parameters of the Hamiltonian (^ (e.g., t±, Vi) may 
be explained by the competition between the electron 
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hopping and the electron-electron interactions, which re- 
duces the gap in the spectrum of isospin excitations. If, 
again, we would leave only the first two terms in the 
Hamiltonian (||), the gap equals \2t± ~ Vi\. Although 
these excitations do not carry charge, they can be ex- 
cited by an electric field applied perpendicularly to lad- 
ders (i.e., in a-direction), which may explain why the 
continuum of relatively low-energy excitations, observed 
in the infrared absorption experiments, disappears when 
the orientation of the electric field is changed from a- 
direction to 6-direction. . The repulsion between elec- 
trons from neighboring ladders (the third term in Eq.(|^)) 
frustrates the zigzag structure, which results in addi- 
tional reduction of Tc- 

Next we consider the spin excitations. From Eq.(^ the 
effective spin-exchange constant in the chain direction, 
Jii, is 



Jn 



/rp rpi rpZrpZ \ 



\\_ IrpXrpX 

^ K-'-n-'-n+l' 



rpyrpy \ 



(4) 



where (. . .) denotes the thermal and quantum average. 

From Eq.(^ it is obvious that CO would decrease the 
effective spin-exchange constant J|| . In particular, if the 
electrons in neighboring rungs n and n -I- 1 of a ladder are 
completely uncorrelated, {T^T^,+i + TyT^^-^) = i, while 
for the extreme zigzag charge ordering ( (2T^) — (—1)"): 
(Tn-Tn+i — T^T^+i) = and the spin-exchange constant 
is zero. In general, the charge order parameter 77 does not 
have the maximal amplitude, so that CO would result 
only in a partial reduction of Jy. According to Eq.(||), 
the spin stiffness in the chain direction becomes smaller 
as temperature decreases, which could be observed in 
neutron scattering experiments. 




FIG. 3. This figure demonstrates that the zigzag CO makes 
the spin-exchange in neighboring pairs of electrons {AB and BC) 
inequivalent. The inequivalence is caused by the difli'erence in the 
occupation of V sites in neighboring ladders. 

More importantly, the zigzag charge ordering opens 



a spin gap, which in this picture is a secondary effect. 
There are, in principle, two scenarios for the spin gap 
opening. First, as one can see from Fig. ||a, the CO 
results in the appearance of pairs of nn vanadiums cou- 
pled by the hopping amplitude t^y, [cf. Fig. For 
large txy these interladder pairs would form spin singlets, 
which could explain the spin gap in NaV205 How- 
ever, we think that this is rather unlikely. On the one 
hand, as we mentioned above, the band structure cal- 
culations give rather small value of t^y ~ 0.012eV [1). 
On the other hand, even with larger t^y (as, e.g., in 
Rcf. [p^) there are several competing mechanisms of nn 
exchange: an antiferromagnetic one, due to direct d-d 
overlap and an exchange via 90° V-O-V path, which ac- 
cording to Goodenough-Kanamori- Anderson rule is fer- 
romagnetic. Even the sign of the resulting interaction is 
unclear: Horsh and Mack obtained a small antiferromag- 
netic nn exchange while the LDA-fU calculations 
for the structurally similar material CaV205 give a fer- 
romagnetic nn interaction Furthermore, in NaV205 
there is a third mechanism of exchange resulting from the 
circular motion of electrons over V-V-V triangles (e.g., 
along the CDE triangle in Fig. ||). The sign of the lat- 
ter exchange coincides with the sign of the product of 
three amplitudes of electron hopping along the sides of 
a triangle [cf. [|9|]. As the later sign is negative the 
corresponding spin-exchange is ferromagnetic. 

Thus, whether the interladder singlets are formed or 
not, is still not clear. It is, however, obvious that the 
zigzag CO results in an alternation of the exchange con- 
stants along the spin chains formed by vanadium lad- 
ders. Indeed, the exchange interaction between the ions 
A and B (see Fig. ||) goes via a plaquette having V^^ 
both above and below it, whereas for B and C ions the 
corresponding positions are occupied by V'^^ ions. This 
inequivalence in the occupation of V sites in neighboring 
ladders produces a difference in the shifts of the corre- 
sponding oxygen and vanadium ions, thereby affecting 
the electron hopping amplitudes and, ultimately, result- 
ing in the alternation of the exchange constants along the 
ladder direction. 



J||(n,n+1) = J„(l-f (-1)"^), 



(5) 



which opens the spin gap 



■■. The detailed calcu- 



lation of the alternation amplitude 5 induced by the CO 
is a formidable problem requiring the knowledge of the 
shifts of oxygen and vanadium ions, the modification of 
electron transfer energies, etc. One can, however, easily 
see, that (5 is a linear function of the charge order param- 
eter 77 — (— )"(2r,j), so the temperture behavior of the 
spin gap is quite similar to that in SP systems, where 5 
is proportional to the amplitude of the lattice distortion 
(unless, of course, the CO transition is of the first order). 

Though the opening of the spin gap is not the main 
driving force of the transition, it also contributes to the 
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stabilization of the charge ordered phase. It may, in par- 
ticular, be important in determining the relative phases 
of charge order parameter ry in different ladders. If, for 
instance, for the phase of the charge order in the bot- 
tom ladder in Fig. || would be opposite, there would no 
doubling of the periodicity in a-direction and no alter- 
nation in the spin-exchange constants along the ladders, 
and consequently, no spin gap. The energy gain due to 
spin-gap opening, by itself, can make the structure of 
Fig. ^ more stable than the zigzag structure without 
the dimerization along a-axis. 

Summarizing, we have shown that the very chemical 
composition of NaV205 - the presence of one electron 
per two vanadium sites - rather naturally leads to a 
charge ordering at low temperatures. We argued, that 
in the ordered phase the occupied vanadum sites form 
zigzags. This CO is driven by the direct Coulomb repul- 
sion and, most importantly, by the interaction with the 
lattice distortions. The CO results in the doubling of the 
lattice period both in a and b directions in agreement 
with experiments. Another consequence of the CO is the 
alternation of the exchange interaction along 6-direction 
(along vanadium ladders), which openins a spin gap. Al- 
though the coupling between spin and charge degrees of 
freedom may also contriubte to the charge ordering, the 
spin-gap opening is a consequence of the charge ordering. 
This picture is consistent with experimental observations 
for NaVaOg. 

We are very grateful to A. Damascelli, R. K. Kremer, 
D. van der Marel, T. M. Rice, W. Schnelle, A. Smirnov, 
A. N. Vasil'ev, and especially to B. Biichner for useful 
discussions. This work was supported by the "Sticht- 
ing voor Fundamenteel Onderzoek der Materie (FOM)" 
and by the European Community through the network 
OXSEN. 



Note added: After completion of this work two preprints 
appeared |^,^ in which similar ideas were put forward. 
Thalmeier and Fulde argued in favor of the CO shown 
in Fig. ||b, after which they still had to invoke the spin- 
lattice coupling to get second transition of a spin-Peierls 
type. Experimentally two close phase transitions were 
observed in some samples of NaVaOs but not in oth- 
ers Furthermore, when two transitions are observed, 
their temperature shifts, essentially, in the same way with 
magnetic field 1^. Thus, it may well be that these multi- 
ple transitions originate from the inhomogeneity of some 
samples, although this question definitely deserves fur- 
ther study. However, as we argued above, the zigzag 
structure of Fig. ||a seems to us much more favorable. 
Besides, in the model of Ref. one completely ignores 
the largest electron hopping amplitude tj_, which would 
result in the random exchange between electron spins 
above Tc, instead of the quasi-one-dimensional behavior. 
The results of Seo and Fukuyama EG] are closer to 



ours; they also concluded that the zigzag structure is 
more favorable. Their arguments, however, are based on 
nn and nnn Coulomb terms only, which is definitely in- 
sufficient for the Coulomb interaction (as we mentioned 
above, the long range Coulomb interaction seems to favor 
the V^^-V^^ chain structure). Furthermore, their sin- 
glets are formed by nearest-neighbor pairs between lad- 
ders, which, as we argued above, is rather questionable. 
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